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. The turbulent Prandtl number has been calculated in the two-loop approximation of the e ex- 

' pansion of the stochastic theory of turbulence. The strikingly small value obtained for the two-loop 

(N ■ correction explains the good agreement of the earlier one-loop result with the experiment. This 

^ . situation is drastically different from other available nontrivial two-loop results, which exhibit cor- 

^ ' rections of the magnitude of the one-loop term. The reason is traced to the mutual cancellation of 

additional divergences appearing in two dimensions which have had a major effect on the results of 
previous calculations of other quantities. 
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I. INTRODUCTION 



>>: 
V 

pH , The method of renormalization group (RG) in the theory of developed turbulence is currently the most developed 
. ■ technical means allowing for reorganization of the straightforward perturbation theory, whose huge expansion param- 
^ ' eter at large Reynolds numbers renders it practically useless. At the same time the physical value of the artificial 
• ^ . expansion parameter e introduced in the RG approach is not small either. For some important physical quantities, 
such as the critical dimension of the velocity and effective viscosity, it is possible to prove with the use of of Galilei 
(-H invariance of the theory that the corresponding series in e terminate at the linear terms. Therefore, for such quantities 
O ^ the RG approach yields exact answers coinciding with the prediction of the phenomcnological theory of Kolmogorov. 
'—'I For other interesting quantities, such as the Kolmogorov constant, skewness factor, turbulent Prandtl number and the 
^ ' like, the series in e, however, do not terminate. In this context, it has been often suggested that with the aid of the e 
^ ] expansions it is not possible to obtain a sufficiently good estimates of numerical values of these quantities, although 
I - until recently ~ there were no calculations extending beyond the first order of the perturbation theory (one-loop 
T-H " approximation) . The two- loop calculation of the Kolmogorov constant and the skewness factor in the inertial range 
, carried out in Ref. confirmed this pessimistic point of view on the whole: the two-loop contribution turned out 
' to be practically equal to the one-loop contribution, although the trend of change of the quantities calculated was 
' ^ , correct, i.e. towards the experimental value from the one-loop result. 

■ In Ref. calculations were carried out for space dimensions d different from d = 3 as well. It turned out that 
■"".^ the relative magnitude of the two-loop contribution decreases with the growth of d, and in the limit d — > cxd is of the 
^ order 10 % only. At the same time in the limit d — > 2 this contribution grows without limit. Such a behavior of the 
coefficients of the e expansion may be related to that their singularities as functions d lie in the region d < 2. The 
nearest singularity at d = 2 is connected with the divergence of some graphs in the limit d ^ 2, which leads to the 
(— I ' appearance of poles in d — 2 in the coefficients of the e expansion, and it is just these graphs which turn out to be 
O I responsible for the large value of the two- loop contribution at d = 3. This feature gave rise to the hope that summation 
of these singularities may lead to quantitative improvement of the results of the e expansion in the real dimension 
. !^ d — ?>. Such a summation was successfully carried out in the framework of the RG method with the aid of the account 
of the additional UV renormalization of the theory in the vicinity of d = 2 . In the resulting " improved e expansion" 
\—{ the low-order terms are calculated in the usual way at d = 3, while the high-order terms are approximately summed 
with the account of their leading singularities in d — 2 (one- loop approximation) , then next-to- leading singularities 
(two-loop approximation) etc. Calculation of the Kolmogorov constant and skewness factor according to this program 
has demonstrated an essential decrease of the relative impact of the two-loop contribution and led to a fairly good 
agreement with the experiment 

In the present paper we shall analyze to what extent the singularities of the e expansion show for another important 
characteristic quantity of turbulent systems, the turbulent Prandtl number. It was calculated in the framework of the 
RG and the e expansion in Refs. 0j 3 with rather good agreement with the experiment 0, 1^0. We have carried 
out a two-loop calculation of the Prandtl number in order to check, whether this agreement is partially coincidental. 

Let us remind that the Prandtl number is the dimensionless ratio of the coefficient of kinematic viscosity vq to the 
coefficient of thermal diffusivity kq. (In the formally identical problem of turbulent diffusion the ratio of the coefficients 
of kinematic viscosity and diffusion is called Schmidt number). For systems with strongly developed turbulence the 
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process of homogenization of the temperature is strongly accelerated, which is reflected in the value of the effective 
or turbulent coefficient of thermal diffusivity. The ratio of the coefficient of turbulent viscosity and the coefficient of 
turbulent thermal diffusivity is the turbulent Prandtl number. Contrary to its molecular analog the turbulent Prandtl 
number is universal, i.e. does not depend on individual properties of the fluid. For the accurate determination of the 
turbulent Prandtl number a set of conditions is required, especially when calculations are carried out in the two-loop 
approximation. Therefore, apart from the formulation of the stochastic problem we shall pay the proper attention to 
this problem as well. 

The present paper is organized as follows. In Sec. ^ we remind the main features of the description of passive 
advection of a scalar quantity in the stochastic theory of turbulence with special emphasis on the careful definition 
of the turbulent Prandtl number within the model considered. Sec. Illll is devoted to the analysis of renormalization 
and renormalization-group equations of the model. In Sec. IIVI details of the two-loop calculation are presented. Sec. 
contains analysis of the results and concluding remarks. 



II. DESCRIPTION OF THE MODEL 



Turbulent mixing of a passive scalar quantity is described by the equation 

dti^ + {^,d,)^ = + f . (1) 

The field ^{'>^, t) in Eq. ^ may have the meaning of both the non-uniform temperature (kq being the coefficient of 
thermal diffusivity) and concentration of the particles of the admixture (in this case is replaced by the coefficient of 
diffusion). The field /(x, t) is the source of the passive scalar field. In the stochastic model of turbulence the field of 
turbulent eddies of the velocity of the incompressible fluid t) satisfies the Navier-Stokes equation with a random 
force: 

dti^i + {(pjdj)(pi = t'oAc^Si - diP + Fi , (2) 

where P(i,x) and Fi(t,x) are, respectively, the pressure and the transverse external random force per unit mass. For 
F Gaussian distribution with zero mean and the correlation function 

{F^{t,x)FJ{t\x')) =5{t-t'){2TT)-'^ j dkPy(k)di.(fc)exp[ik(x-x')] (3) 

is assumed. Here, Pij(k) 5^ — k^kj/k'^ is the transverse projection operator, dp(k) a function of fc = |k| and 
parameters of the model, and d the dimension of the coordinate space x. 

The stochastic problem ([Ql - (O is equivalent to the quantum-field model with the doubled number of fields 
(j) = ip, ip\ ^/;'} and the action 

5($) = ^'Bf ip'/2 + if'i-dtip + i^o^ip - {ipd)ip] + V/[-5t^ + KoA^ - {^d)^ + /], (4) 

in which _Dj? is the correlation function of the random force ^ and the necessary integrations over {i,x} and 
summations over vector indices are implied. In model only correlation functions of the admixture field of the 

form 

(V(xi,ti),i/;(x2,t2)...^/^(x„,t„)V^'(x;,i'i),^'(x'2,i^)...^A'(x;,C)), 

with the meaning of multiple response functions are nonvanishing. The simplest of them is determined by the following 
variational derivative with respect to the source / in Eq. 



G(X - X', i - t') EE (^(X, i)^'(x', t')) I = ^^^^""^ 



i/=o <5/(x',t') 



(5) 



The non-random source field / of the passive scalar has been introduced in action solely to remind of relation (0 
and its generalizations and will therefore further be omitted. 

Model gives rise to the standard diagrammatic technique with the following nonvanishing bare propagators 

{t = ti-t2) 



{ip{hMt2)h - 7^ exp {-i^ok'\t\) , (6) 



2vok^ 

{ip{h)ip'{t2))Q = e{t) exp {-vok^t) , (7) 
mi)^p'{t2))o = 0{t) exp (-^iokh) , (8) 
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in the {t, k) representation. The common factor Pij(k) has been omitted in expressions © - (U)) for simpUcity. 
Interaction in action corresponds to the three-point vertices —ip'{ipd)ip — ip'^Vijs^Pjfs/'^ with the vertex factor 
Vijs = i{kj6is + kgSij), and —^'{ipd)ip = ikjtjj'ipj^p, where k is the wave vector of fields ip' and ijj' . 

Turbulent processes lead to significantly faster than in relations Q) and JSJ attenuation in time of the response 
functions {f(p') and {ipip') due to the effective replacement of the molecular coefficients of viscosity and thermal 
diffusivity by their turbulent analogs. At the same time, however, the simple exponential time-dependence is changed 
as well (and in a different manner for {(pip') and (ipip')), therefore it is necessary to choose a definite way of fixing 
the ratio of the turbulent transport coefficients, i.e. the Prandtl number (or Schmidt number). Henceforth, we shall 
use the following definition. Consider the Dyson equations for the response functions in the wave-vector-frequency 
representation: 

G~l,{k,uj) = r^<p/(fc,cj) = -luj + vok^ ~ E^'^(/c, w) , (9) 
G^]p, {k,uj) = r^^' (fc, uj) = -\Lu + kq/c^ - S^'^ {k,uj) , (10) 

where S are the corresponding self-energy operators, and introduce the inverse effective Prandtl number Me// by the 
relation 

_ r0^/(fc,cj = o) 

"^^^ = r,,,(fc,c. = o) ■ (^^^ 

Further, we shall be interested in the inertial range <^k <^ A. (here, L is the external scale of turbulence and 
the characteristic length of the dissipating eddies) in which the quantity Me// is independent of the wave number k. 



III. RENORMALIZATION OF THE MODEL AND THE RG REPRESENTATION 

The self-energy operators ^^{k,uj) and S^'^(/c,a;) appearing in Eqs. ^ and (|10|l may be found in model (0)) 
in perturbation theory. However, the expansion parameter turns out to be very large for developed turbulence (for 
Ai ^1). The renormalization-group method allows to carry out a resummation in the straightforward perturbation 
theory. To apply it, it is necessary to use in relation (|3Jl "the pumping function" dp{k) of a special form 

dF{k)^Dok^~'^'^'. (12) 

In the infrared region the power function (|12|l is assumed to be cut off at wave numbers k < m = L~^. The quantity 
e > in Eq. p2(l is the formal small expansion parameter in the RG approach with the value e = 2 corresponding to 
the physical model. 

The usual perturbation theory is a series in powers of the charge go = Dq / Vq dimensionless at e = (logarithmic 
theory). At e — > ultraviolet divergences are brought about in the graphs of the perturbation theory which show 
in the form of poles in e. Due to Galilei invariance of the model divergences at c? > 2 are present only in the one- 
irreducible functions {(p(p') and {ipip') and are of the form (p'Aip and il^'Aip. At d = 2 also the one-irreducible function 
{f'if') diverges. For d > 2 the renormalized action may be written as 

Sr^ = i ^'Df^' + ^'[-dt^ + j/Z„A^ - {^d)ip\ + ^''[''dt^ + uiyZiA^ - iipd)i>] . 

It is obtained from action Q by the multiplicative renormalization of the parameters of the model: 

i^o = i^Z^, go ^ gi-i^^Zg, uq ^ uZ^, Zu = ZiZZ^, Zg = Z;;^ (13) 

with two independent renormalization constants Zi, and Zi . The quantities v and g in Eq. (|13|l are the renormalized 
analogs of the coefficient of viscosity and the coupling constant (the charge g being dimensionless). The renormalization 
mass fi is an arbitrary parameter of the renormalized theory, and the pumping function dp{k) 112|) determining the 
correlation function of the random force Dp ^ is assumed to be expressed in terms of the renormalized parameters: 

dp{k) = g^vlk^-"-^^ = g^,^^v^k'-''-^^. 

The dissipative wave number A is determined by g^ according to the relation A = gQ^^'^ . It may be also estimated by 
the quantity /i. Thus, the inertial range we are interested in corresponds to the condition s = fc/^ <C 1. 
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In the scheme of minimal subtractions (MS) used in the following the renormalization constants have the form of 
the Laurent expansion 1 + poles in e 

oo oo n 

Z = 1 + ^ afe(g, u)e-'' - 1 + E E • (14) 

k=l n=l k=l 

For Zi, at (i = 3 in Ref. '3| the following expression was obtained 

Z. = l + ^ + 0{g), a,,=-^^^, 5. = ^, (15) 

where Sd = 271^^/^ /T{d/ 2) is the area of the d-dimcnsional sphere of unit radius. 

The correlation functions of the renormalizcd theory do not contain poles in s. This feature, however, does not 
solve the problem of finding the infrared asymptotics s = — > 0, because the corresponding perturbation theory 
is a series in the parameter s^^^ growing without limit in the region we are interested in. The problem is solved by 
passing to the RG representation. To use it for the response functions © and Hl()(l , rewrite them in the renormalizcd 
variables in the form 

r^^,{k,LU = 0) = vPR^{s,g), 

T^^,{k, Lj ^0) = uvk^R^{s,g,u), (16) 



where the dimensionless functions R^p and R^ of dimensionless arguments s, 5, and u are given by the expressions 

R^{s,g) = Z^- 



Yi^,^{k,uj = 0) 



R^{s,g,u) = Z,- -■ (17) 

The RG representation for functions (|16|l is determined by the relations 

T^^>{k,ui 0) = vk'^R^p{s l,g) , 

r^^' (fc, w = 0) = uDk^R^{s = l,g, u) , (18) 

where g = g{s, g), D = v{s, g, v), and u — u{s, g, u) are invariant variables satisfying RG equations of the form 

[ - s(9s + Pgdg + Pudu - Iv vd^] b{s, 9, u) = 0, 

and normalized by the conditions g(l, g) = g, v{l, g, v) = v, and u(l, g, u) = u. The RG functions (3 and 7 are defined 
by the renormalization constants according to the relations 

Pg{g) = ^l^^\^g = g{-2e + Z-i^), /3u{g,u) = nd^\^u = ^(71 -71,) , 

71,(5) = I pin , 7i(g,u) = ^a^lglnZi, (19) 

where M^^Iq denotes the operator /ic?^ acting at fixed bare parameters go, vq, and uq. The last equalities for the /3 
functions in Eq. p9fl are a consequence of the connections between the renormalization constants in Eq. (|13|l . 

As shown in the one-loop approximation in Ref. 0,00, the invariant charges g{s,g) and u{s,g,u) in the limit 
s — > tend to the infrared-stable fixed point: g{s,g) g* = 0(e), u{s,g,u) ^ = 0{£^), and the invariant 
viscosity has the powerlike asymptotic behavior 

Dok-^^V'^ fDok-^^^'/^ 



Thus, the expression for the effective inverse Prandtl number l|ll(l in the inertial range predicted by the RG represen- 
tation with the account of relations H16|l and (|18|1 is 

i?^(g = 1, g*: M*) 

Ueff = p . _ -I „ ^ • (20) 

H^[s — 1, g*j 
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IV. TWO-LOOP CALCULATION OF THE PRANDTL NUMBER 

The expansion of the functions i?^ and (|17|l in the couphng constant g is of the form 



Rip — 1 



Roi 



s 



A^{u)s' 



-2e 



(21) 



Here, the quantities A^p and are determined by the one-loop contribution to S,^'^ and S^'^, whereas the coefficients 

a"ii and a^^i of representation (|14() of the renormahzation constants Z^, and Zi are found from the condition of UV 
finiteness of expressions (|21|l . Substituting relations H21|) in Eq. H2U|) we obtain 



e// = "* {l + ["v ~ aV'('"*)]3* + O {gl)} 



(22) 
(23) 



Bearing in mind that — 0{e), we see that to find Ueff at the leading order of the e expansion, it is enough to know 
the charge u* in the one- loop approximation. At the second order, apart from the more accurate values of and g*, 
it is necessary to calculate the coefficients and a^(u*) of the expansion of the scaling functions H17|l and (|21|l at 
the leading order in e as well. 

The location of the fixed point (g*, u^) is determined by the conditions Pg{g*) — /3«(<7*,u*) = 0. The nontrivial 



fixed point with 5^ ^ is infrared stable and from Eq. (|19|1 the relations 



7i(g*,w*) 



2e 

y ' 

2e 

y 



(24) 
(25) 



follow at this fixed point. 

The UV-finiteness of the RG functions 7(3, m) from Eq. (|19|l allows to express them in terms of the coefficient of 
the first-order pole in s in expression H14|l for the renormahzation constants: 



7 = iPgdg + Pudu) InZ = ~2gdgai . 



(26) 



The renormahzation constant Z^, at the second order of perturbation theory and the corresponding expression for 7^ 
were obtained in Ref. For g^ from Eq. (|24|1 the result is: 



(27r)° 



(27) 



where 



-1.101, 



3(d- 2) 



+ c + Oid-2), d 



From previous analyses the renormahzation constant Zi is known in the one- loop approximation only Q: 



Zi = 1 



flu .9 



4 + ^) {gS,r + 0{g^) 



,(1) 



id - l)Sd 
4du(l + u) ' 



(28) 



(29) 



Calculation of the contributions A{u) and B{u) of order g^ is presented below [it should be noted that, like the 
one-loop factor a\^i , the two-loop coefficients of the poles in e in representation (|29|l are non-polynomial functions of 
u] . According to Eq. , the RG function 71 corresponding to Eq. is 



71 



id ~ l)gSd 
2dw(l + u) 



^gSdf + 0(5' 
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FIG. 1: The one-loop self-energy graph for Ti^i^. The lines correspond to propagators © and ||7J. Slashes denote the 
end carrying arguments of the field (p'; plain end carries the arguments of the tp field. Vertices correspond to the factor 




FIG. 2: The one-loop self-energy graph for E^/^. The lines correspond to propagators Q and Q. Slashes denote the end 
carrying arguments of the field plain end carries the arguments of the field ip or tp. Vertices correspond to the factor 



iki 



Iterative solution of Eq. (|25|l with respect to u with the account of relation H27|) yields 



,(0) 



,(0) 



2{d + 2) 



,(1) 



2{d + 2) 



1 



128(d + 2)2 / 



(30) 
(31) 



Substituting relations H30() and (|31|l in Eq. (|22|l and taking into account Eq. (|27|l we obtain 



,(0) 



1 



,(0) 



l + 2u 



(0) 



128(<i+2) 
3(d-l)2 



,(0) 



8(d-f 2) 
3(d-l)5d 




(32) 



We now turn to the calculation of the constants B, a.^ and which determine the Prandtl number. In the one-loop 
approximation the quantities Y^^'ip and S^'^ are represented by the graphs depicted in Figs. ^andlJl respectively. 

In the one- loop self-energy graphs of Figs. ^ and 13 the lines correspond to propagators Q and © with the 
convention that ends with slashes corresponds to arguments of the fields ip' and ip' , plain ends of Lp and tp. Vertices 
in Figs, n and 121 correspond to the factors Vijs = i{kjSis + kgSij) and ikj, respectively. Upon contraction of indices, 
integration over time and introduction of dimensionless wave vector (in units of the external wave vector p) in the 
integrals we obtain 



2{d- 1) J {2Tr)d 



dk /c2-<i-2£(^ „ C2)[2fc3^ _ _ 3)^2 _ 2k{d - l)e -{d-1) 



{2k^ + 2k^ + l){k'^ + 2k^ + l) 



dk 



{2ttY {l + u)k^ + 2uk^ + u 



(kp) 



(33) 
(34) 



The integrals H33I) and H34|) are UV divergent in the limit £ ^ 0, the residue at the pole is readily found by selecting the 
asymptotic at fc — > oo contributions to the integrands and discarding the inessential region of integration fc < 1. Thus, 
for the coefficients and together with the renormalization constants and Zi chosen to cancel divergences in 
expressions (|23|) we find 



1 



Aid - l){2TrY J k 
1 



dk 



-1 

2u(l-Ku)(27r)° 



dkil~e)m~d + 3~6e), 

'dkii^e), 



dk 



k 

I' 
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Replacing the integral over directions of the unit vector k by the average over its directions J dk... = Sd{---) and 
taking into account that 

we arrive at result (|15|l for a'^ii and (|29|l for a^\\ In view of the preceding argumentation, the coefficients a,p and 
in Eq. H23|) at the leading order in e may be written as 

oo ✓ 

1 f^, f^ir, ^2nJ 2fc[2fc3^-(d-3)/c2-2fc(d-l)e-(d-l)] 



4(d- l)(27r) 



^Jdkjdkd U<j (2fc2 + 2fce+l)(P + 2fce + l) 

6'(A:- l)(2/fcC-«^ + 3-6^2) 



, (36) 



" 2ii(27r)'i 



At d = 3 from relations H36|) and H37|l we obtain 

_ oo 1 

S: 



k 9{k-l) 
(1 + ?i)fc2 + 2Mfc^ + u k{l + u) 



:^ jdk jd^i-e) 



-^3 



2k[k^^ - 2k^-l] 



dk /d^(l-C^) 



(2fc2 + 2fc^+l)(/c2 + 2fc^+l) 

k 9{k-l) 



-9ik-l){^-'4. 



{1 + u)k'^ + 2uk^ + u k{l + u) 



d ~ 3, S3 — —T , 
27r^ 



(37) 



(38) 
(39) 



Numerical evaluation of integrals (|38|l and H39() with u — ul"'' from Eq. (|30|l yields 

= -0.047718 S3, = -0.04139 ^3 . 



(40) 



It is convenient to find the two-loop contributions to the renormalization constant Zi from the condition that the 
quantity from Eq. (|17|l is UV finite in the limit fc ^ 0. In terms of the reduced quantity 



lim 



T,^',^{uj = 0, k) 
uvk'^ 



this condition may be cast in the form 



Zi(£)-S(e) = O(£0). 



(41) 



(42) 



The limit fc in expression (|41|l does exist, provided the IR regularization of the graphs has been taken care of. 
In the MS scheme renormalization constants do not depend on the method of such regularization. With our choice of 
the pumping function (fT^ it is accomplished by the cutoff of the propagator {v5</')o © at fc < m. 

Let us choose further the wave vector of integration such that in the lines {f(p)o it flows alone (for the graphs S^',^ 
such a choice is always possible). Then integration over all the wave numbers will be carried out within the limits 
from m to 00. 

The one-loop contribution to S is determined by the graph of Fig. El as: 



S(i) 



2e 



2u{Z^ + uZi)Z^ J {2Try- 



dk (1 ~ ^^)9{k ^ m,) 



5M 



2s 



dk 



2u(Z^ + uZi)Z^{2nY J k<^+^^ 



dk{\~e)^- 



2e 



dk 



2w(Z^ + mZi)Z^ J k<^+^^ 




2:,= - 




L = ^ 

^5 2 





^7= -. 





f^--f ---+^— 



FIG. 3: The two-loop self-energy graphs for E^,/^. The lines correspond to propagators 0, Q and Slashes denote the 
end carrying arguments of the field tp' or plain end carries the arguments of the field (fi or tp. Vertices correspond to the 
factor Vijs = i^kjSis + ksSij) or Vijs = ikj. 



which, together with relations H15|) . H29|) . and H35|) yields 

4:eu{Z^ + uZi)Zt,d 4eit(l + M)d 

where t = m/fj,. Extracting the pole contributions in e from expressions (|43ll we obtain 



ua[\^ + (2 + u)aii 



gSd 

e{l + u) 



+ 0{g'), (43) 



Ed) 



4£u (1 + u) d 



taiV + (2 + ")an 



gSd 

(l + u) 



21nr 



0(e°) 



(44) 



Substituting relation H44|) in Eq. H42|) and requiring cancellation of pole contributions in the linear in g approximation, 

we return to expression H29|l for aj^-*. The terms of order are required for the calculation of the renormalization 
constant in the two-loop approximation. 

The two-loop contribution E^^) to the self-energy operator S^/^ is determined by the sum of the graphs depicted in 
Fig. 13 [normalization according to Eq. (|1T|| is imphed]. When substituting propagators ©, Q) and - expressed 
in terms of the renormalized variables - in the graphs of Fig. |21 it is possible to put = Z\ = \ with the necessary 
accuracy. Contracting indices and integrating over time we obtain 



(.9^d)> 



2,,4e 



dk 



dq 



(1 - e) Jn 



1, 2. 



(45) 



where 



and 



Ji = 2q' 



J2 = —zkq , 



(46) 



2,,4e 



96uv 



dk 



^l+2e 



k^ + 2kq^ + g2 



n = 3...l 



(47) 
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with 



J3 = fc(fc3 + 2k^q£_ - q\) 
(fc3 + 2Pq^ - q^^) 



J5 
J7 



{k'^ + kq^ + q^)ivk^ + kq£_ + q"^) vk^{vk'^ + kq^ + q"^) k^ {k^ + kq^ + q"^) 



k{k^ + kq^ + q"^) ' 

fc2[fc^ + g-* + fcge(fc^ + g^)] 



2kq£_ 



f (fc2 + -|- ufcg^ 



+ 



A;2 + fcg^ + g2 



kq^jk^ - q') 
2{k'^ + kq^ + vq"^) \vq'^ 

P{2k^ + 3k^q^~q^^) j 1 



(48) 
(49) 
(50) 
(51) 



vk^{vk^ + kq^ + q^) {k^ + kq^ + q'^){vk^ + kq^ + q^) 



k{2k^ + 3fc2(7e - q^O 



2(fc2 + kqC + g2)[„(/c2 ^ ^2) ^_ ^jfcq^] 
Integrals - may be represented as 

E,(e) = M 



4e 



fcl+2? 



dq 



de/.(e,fc/9), 



(52) 
(53) 

(54) 



or, after the corresponding stretching of integration variables, as 

1 



Me) 



dk 



dq 



ljl+2e J ql+2e 



dih{Lk/q). 



or, finally, as 



-4e 



Me) 



dk 

/fcl+2' 



dU^itk/q). 



We are interested in the coefficients of the pole contributions to Si(e): 

^. = ^ + ^ + 0(e°), I].(e) = ^+ ^'-'"'^"% 0(s°). 



(55) 



(56) 



(57) 



For the functions fi{z^ k/q) with i — 2, 5. ..8 the equations fi{z, 0) = fi{z, 00) = hold revealing that integrals over k 
and q in Eq. H54|l are separately convergent, so that the divergence at £ —> in the corresponding Ei is brought about 
by the region, in which k and q tend to infinity simultaneously. As a consequence, the second-order pole is absent in 
such E.;: Gi = for i = 2, 5. ..8. 

For E,j with i ~ 1,3,4 fi(z, 00) = as before, which means absence of divergence in the integral over k in Eq. 
For these graphs, however, fi{z,0) = const 7^ 0, so that the integral over q diverges at £ — > leading to the 
appearance of the pole of second order in the full integral. 

Expressions H55|l may be simplified with the use of the identity 



E, = 



4£ 



(58) 



following from Eq. H5t)|) . Calculating the right-hand side of Eq. (|58|1 with the aid of relations (|55|1 and introducing 
the dimensionless integration variables, we obtain 



Me) - 



4£ 



dn 



(59) 
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TABLE L Residues of the first-order poles in e of the dimensionless integrals I56II corresponding to the two-loop graphs of 
Fig. El 



i 


1 


2 


3 


4 


5 


6 


7 


8 


bi ■ 10^ 


0.1099 


0.0944 


0.8691 


0.0057 


-3.9382 


0.0672 


-1.9647 


0.5899 



This operation has reduced the number of iterated integrations and allowed for explicit extraction one pole in e. For 
i = 2, 5. ..8 the integral in Eq. is finite for e = and determines the residue of the first-order pole: 

C30 1 

a, = 0, b, = -^j^ jd^[U{i,n) + f,{^Mn)], z = 2, 5...8. (60) 
1 -1 

For Si(e) with i = 1,3,4 the coefficient of the second-order pole is obtained by the replacement of the function 
+ 1/'*)] the integrand in Eq. (|^ by its hmiting value at k ^ oo: cjo) -|- 0) = 0) [we 
remind that /^(z, oo) = 0]. Then integration over n becomes trivial, which yields 

1 

a. = l Jdtm,0), z = 1,3,4. (61) 

The remaining integral with the change fi{£,, k) — > [fiH, «) ~ fiiS,, 0)] is finite at e = and determines the residue of 
the first-order pole: 

oo 1 

I— I dz[U{z,n) + fdzM^)- h{zM, » = 1,3,4. (62) 



4 ./ K 

-1 



Let us write condition H42() at order for d ~ i. With the use of the corresponding terms of the one- loop contribution 
()44|l . the summed two-loop contributions (|57|l and expression H29(l for the renormalization constant Zi, we obtain 



A B 1 



e Qeu{l + uy 



8 



.a« + (2 + H)4r'J (^-21nr)+^(g + ^'~';'^"M . (63) 



With the aid of expressions and if^H) in Eq. it is not difficult to find 

h) 



72u(l + u)3' ''^ 480m(1 + u)2' 480^(1 + u) ' 

Substituting these values in Eq. 163|) and taking into account relations H29|) and l|15|) for a^^-* and a^^^ , we see that the 
terms with Inr in Eq. I|63|) are automatically cancelled (as a consequence of renormalizability of the model), whereas 
for the coefficient A of the second-order pole we obtain 

-I- 9u -I- 16 



A^ 



720w(H-u)3 



For the coefficients bi numerical integration of expressions ijnni) and (inU with u — ui from Eq. (|30|l yields the 
results quoted in Tabled which for the coefficient B in Eq. H63|l lead to the value 



5 

B (uf^^ = '^bi = -4.1666 • 10-3 . 



Substituting this value in Eq. (|32|l as well as and from Eq. H40I) and A from Eq. H28|) . we obtain the final 
expression for the effective inverse Prandtl number: 



ueff = u'^°\l-0m58e) + O{e^), = ^/^^^ — 1~ 1.3930, d = 3. (64) 

At the physical value e = 2 this yields for the turbulent Prandtl number Pr^ the result: 

Pr|°^ ~ 0.7179 , Prf - 0.7693 (65) 
in one-loop and two-loop accuracy, respectively. 
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V. 



CONCLUSION 



The main conclusion to be drawn from the two- loop value of the effective inverse Prandtl number (|64|l obtained in 
the present paper is that the correction term is strikingly small. Even at the real value e = 2 it is only 7% of the 
leading contribution. Apparently this is the reason of the favorable comparison of the one-loop value of the turbulent 
Prandtl number 0.72 with the experiment: the recent circular-jet result 0.81 ± 0.05 |^ is corroborated by the 

experimentally recommended value 0.8 for modeling and once more confirmed the range 0.7 - 0.9 of measured values 
[jij which has been put forward, however, quite a while ago In view of these numbers it may be concluded that the 
already fairly good one-loop result is improved by the two-loop correction, whose account leads to the value 0.77 
for the turbulent Prandtl number. At the same time this result is somewhat unexpected - similar two-loop corrections 
to the Kolmogorov constant and the skewness factor are large ^ . 

In the results obtained there are, however, also significant features common with the calculation jj. From Tabled 



of coefficients bi we see that 65 has the largest value while the value of the whole sum B (ui j = X)i=i close to 



that of 65 alone. The graph S5, which gives rise to this coefficient, is the only one of the two-loop graphs of Fig. |3 
possessing a singularity in space dimension d = 2. Exactly the same situation was encountered also in Ref. Jj in the 
two- loop calculation of the constant A ~ —1.101 in Eq. H28|) . It is a rather unexpected observation that for terms in 
the factor M = A — 8005/3 ~ 0.010 from Eq. H32() {d ~ 3) an almost complete numerical cancellation takes place. 
To clarify the situation, we calculated asymptotics of the graph S5 at d 2 with the result: 65 ~ — l/1024((i — 2). 
Substitution of this expression together with the analogous one for A from Eq. 128|) in the factor M{d) reveals that 
the singular in c? — 2 contributions indeed cancel, so that this factor and with it the whole two-loop contribution to 
the Prandtl number are finite at rf = 2! It should be noted that also in the second term of the two- loop contribution 
to Ueff [see Eq. H32(l ] a significant decrease in magnitude takes place in the difference — of the contributions 
from renormalization of viscosity and thermal diffusivity compared with the magnitudes of these terms separately. 

Thus, our results complement the conclusion made in Ref. 1]. In the two- loop approximation the main contribution 
is due to graphs having a singularity at c? = 2 and it is necessary to sum such graphs. For quantities in which this 
singularity is absent the two-loop contribution is relatively small. 
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